Coulomb drag at v 



Composite fermion pairing fluctuations 
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We consider the Coulomb drag between two two-dimensional electron layers at filling factor v — ~ 
each, using a strong coupling approach within the composite fermion picture. Due to an attractive 
interlayer interaction, composite fermions are expected to form a paired state below a critical tem- 
perature T c . We find that above T c pairing fluctuations make the longitudinal transresistivity p D 
increase with decreasing temperature. The pairing mechanism we study is very sensitive to density 
variations in the two layers, and to an applied current. We discuss possible relation to an experiment 
by Lilly et al. [Phys. Rev. Lett. 80, 1714 (1998)]. 



Double layer systems in the quantum Hall regime 
exhibit a rich variety of physical phenomena, includ- 
ing transitions between compressible and incompressible 
states, as well as interlayer coherence jjj. A powerful 
probe to such systems is a Coulomb drag experiment, in 
which a current I\ is driven through one layer, and the 
voltage V% is measured on the other layer. The transre- 
sistivity is defined as p D = —V2/I1 for a square sample. 
A Coulomb drag measurement probes the interlayer in- 
teraction, the single layer response, and interlayer corre- 
lations. 

We concentrate here on the Coulomb drag between 
identical layers at filling factor v — i each, which has 
been recently measured by Lilly et al. pi . The weak cou- 
pling limit of this problem was considered using a com- 
posite fermion picture |5]Q] and an electronic picture || , 
leading to three main results. First, p D is expected to be 
much larger (typically by 3-4 orders of magnitude) than 
at zero magnetic field, due to the slow decay of density 
fluctuations B. Second, the low temperature dependence 
is p D oc T 4 / 3 . Finally, within a weak coupling the- 
ory p D vanishes at T — 0. The observed low temperature 
behavior || seems to contradict the last conclusion, and 
indicates the need to go beyond weak coupling theory. 

In this paper we consider the Coulomb drag using a 
strong coupling approach within the composite fermion 
picture. Part of the interlayer interaction between 
composite-fermions is attractive Q, and as a result a 
composite fermion paired state is expected to form be- 
low a critical temperature T c 0. In electronic terms, 
such a state is an incompressible quantum Hall state of 
the double layer system, and is expected to have a large 
transresistivity, similar to that discussed in Ref. ||. Ex- 
perimentally this state is not observed, and p D is always 
much smaller than the single layer resistivity p xx || . Mo- 
tivated by this experimental finding, we consider here the 
contribution of fluctuations of the paired state above T c . 
Our central result is that by including the contribution 
of pairing fluctuations the transresistivity becomes 
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where a 



h T — T c 

0.12, and p^ is the weak coupling con 



tribution §. This expression is valid provided that 
T — T c <C T c , but p D <C Pxx- Unlike po \ the first term in 
(JjJ) increases with decreasing temperature. While a sim- 
ilar term affects p xx as well, p D is much more sensitive to 
it, being typically three orders of magnitude smaller than 
Pxx ■ Within the mechanism we study, the drag voltage is 
always parallel to the current, i.e., there is no Hall drag. 

In the composite-fermion approach [S, the system 
of electrons is transformed by a suitable Chern-Simons 
transformation to a system of composite fermions in an 
average zero field, carrying flux lines attached to them. 
In the double layer system, if the interlayer spacing d 
is large, a single-layer transformation is to be performed 
on each layer separately, leading to the Hamiltonian, in 
units where % = c = 1, 

i 
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Here, ipj(r) is the creation operator for a composite 
fermion at point r in layer i, mb is the electron's band 
mass, 6pi(r) = ■0j( r )V , i( r ) — n is the density (with 
the average value n subtracted), <5a,;(r) is the Chern- 
Simons gauge field (with the external vector poten- 
tial subtracted), and the colons stand for normal or- 
der. The density and the gauge field are related by 
An5pi(r)i — V x 5a.i(v). The Coulomb interaction is 
given by %(r) = e 2 /ey / r 2 + d 2 (l — Sij), where e is the 
bulk dielectric constant. Since the Chern-Simons trans- 
formation is carried out separately on each layer, the elec- 
tronic transresistivity (both longitudinal and Hall com- 
ponents) is identical to that of the composite fermions. 
At v — i, composite fermions experience no magnetic 
field on average. Consequently, their Hall transresistiv- 
ity vanishes, and so does the electronic one. We find this 
result to hold also around v = -|. 

The composite fermions interact through the Coulomb 
interaction and through an intralayer Chern-Simons in- 
teraction introduced by the transformation (0). Within 
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the random phase approximation, the most singular ele- 
ment of the gauge propagator is between anti-symmetric 
transverse components. This element mediates an inter- 
layer interaction between transverse currents fj , 
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Here, composite fermions of momenta ki, k 2 are scat- 
tered to ki+Q, k 2 — Q, the energy transfer is f2, the Fermi 
momentum is k F — y/Aim, and the limits Q <C qk F /m* 
and q <C k F are used. The energy scale To is, in the 
limit of strong screening, To = ire 2 nd/e ||. The ratio of 
To and E F ~ e 2 k F /e is proportional to k F d. We assume 
T ^> E F ^> T. We use a dressed interaction vertex, 
inversely proportional to the composite fermion effective 
mass m* . The contribution of all other components of 
the gauge propagator is neglected. The strong singular- 
ity of (||) in the limit Q, — > 0, which results from the 
slow dynamics of charge relaxation in the v — 1/2 state 
p), poses the main difficulty in the present calculation. 

From Eq. (BH it follows that the interlayer interaction 
between fermions moving in opposite directions (k 2 = 
— kx) is attractive p|. The only bare interlayer interac- 
tion is the electrostatic one. Thus, this attractive in- 
teraction may be understood in terms of the electronic 
Coulomb interaction: a composite fermion in momentum 
state k carries a current j oc k. In electronic terms, such 
a current is accompanied by a chemical potential gradi- 
ent (2h/e 2 )z x j. The v = \ state is compressible, so that 
a chemical potential gradient leads to a density gradient, 
which is proportional to z x k. The attraction between 
a composite fermion of momentum k in the first layer 
and a composite fermion of momentum — k in the second 
layer results, then, from the Coulomb attraction between 
the density gradients accompanying the two composite 
fermions, which are oppositely directed. 

The composite fermion attraction is expected to lead to 
interlayer Cooper pairing below T c M. Above T c , pairing 
fluctuations enhance the composite-fermion conductivity 
tensor, and hence p D . As in conventional superconduc- 
tivity JTofl , there are Aslamazov-Larkin |nj and Maki- 
Thompson |L2"| contributions. Pair breaking mechanisms 
are significant in this problem (see M and below). Con- 
sequently, the Aslamazov-Larkin term dominates close to 
T c . Physically, it corresponds to the enhancement of the 
conductivity due to the formation of composite fermion 
Cooper pairs by thermal fluctuations. 

We now outline the calculation leading to Eq. (Q). A 
detailed presentation will appear elsewhere. We use the 
finite temperature Matsubara Green function method. 
To treat the divergences resulting from the singular in- 
teraction (j^) at finite temperatures, we separate the in- 
teraction, depending on the bosonic Matsubara energy 
transfer (fi m = 2nmT), into a static part, W stat , when 
VL m = 0, and a dynamic part, U dyn , when Q m ^ 0. The 



infrared divergences result from the static part, which 
may be viewed as analogous to a singular disorder po- 
tential that is identical in both layers. These divergences 
are exactly canceled in the calculation, a result which is 
a manifestation of Anderson's theorem. 

We use the following notation: Calligraphic letters 
(e.g. U) denote matrices, with indices corresponding to 
both momentum and Matsubara energy. Matrix multi- 
plication is defined with the accompanying factors, 

v-^ f d 2 k 

(A * ff)kro,k'm' = —T y j j (2 \2 ^km,kAm,k'ra' " (^) 

The trace over a matrix of this type (denoted by tr*) and 
a multiplication with a vector, which is marked with an 
arrow, are similarly defined. 
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FIG. 1. Diagrams for (a) The Aslamazov-Larkin contribu- 
tion to the composite-fermion transconductivity (b) The 
particle-particle ladder K, (c) The static ladder C. Here, W stat 
is marked with a dashed line, and U dyn with a dashed wig- 
gly line. The Green function is denoted by a bold line. The 
vertex function qjT in (a) is derived from Eq. (|l8|). 

The diagrams involved in the calculation of the 
Aslamazov-Larkin contribution are given in Fig. @. The 
composite-fermion transconductivity is given by 
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Here, /C(q, ui m ) is the particle-particle ladder [Fig. |l](b) 
and (c)], with total incoming momentum q and total 
incoming (bosonic) Matsubara energy oj m . The vertex 
function is given by c\J , where J , unlike /C, is regular 
in the long- wavelength low-energy limit at T = T c . The 
external frequency is flm- The calculation is carried out 
to leading order in T — T c . 

Because of the separation of the interaction into a 
static and dynamic part, the equation for the particle- 
particle ladder K, is also separated into two parts. The 
integral equation for K. [see Fig. |l](b)] is given by 
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AC(q, u m ) = W dyn (q) + W dyn (q) * C(q, w w ) * /C(q, o; m ) . 

(6) 

where the static ladder C is the solution of the integral 
equation [see Fig. |l](c)] 

C(q, u) m ) = 0(q, w m ) + G(q, uj m ) * W stat (q) * C(q, w m ) . 



(7) 



Here. 

5km,k'm'(q,Wm) 



,<5 2 (k-k') 



(8) 



G(k + q, e m + w A )G(-k, -e m ) . 



where e m = (2m + 1)ttT is a fermionic Matsubara en- 
ergy. The static ladder is analogous to the Cooperon in 
disordered systems. For the Green function, G(fc, e m ) = 
[ie m — £fe + S(/c,e m )] _1 , we use the self-consistent Born 
approximation for the static interaction, 



m cm) =T J = 0)G(fc', e m ) . (9) 



Eq. (g) results from the observation that the most sin- 
gular term in the intralayer interaction is —U\\. The 
factor T appeared from the separation of the Matsubara 
sum. The contribution of U Ayn to the self-energy is non- 
divergent. Thus, its inclusion here is not essential. The 
angular averaging of U atat is, for small \k — k'\, 
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where k and kl are assumed be close to k F , and 9 is 
the angle between k and k'. The divergence due to the 
interaction is thus logarithmic. To regularize the diver- 
gence, the singular term l/\k — k'\ may be replaced with 
tanh[A(fc — k')]/(k — fc'), and the limit A — > oo taken at 
the end of the calculation. In our case the terms involving 
ystat are cance i ec i exactly. 

We first describe the solution of the particle-particle 
ladder [Eqs. §)-(§)] at q — ui m — 0, where the pairing 
instability emerges at the critical temperature T c . The 
static ladder C has a singularity at e = (after analytic 
continuation ie m — > e + iS), where resonant scattering 
occurs. The terms contributing to the singular part C s 
are those with multiple scattering lines. As a result, C s 
is separable in the momentum variables: 
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where N is the density of states at the Fermi energy and 
the energy is measured from the Fermi surface. 

The particle-particle ladder /C(q = uj m = 0) is singu- 
lar at T = T c . To calculate the singular part JC S [using 



Eq. (H)], an angular averaging over the Fermi surface of 
U dyn {(\ = 0) is needed [cf. (pi)]. For m ^ m! it is 
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This result assumes To(|fc — k'\/k F ) 3 <C |£2 m |, an assump- 
tion justified a posteriori by Eq. (|l5|). After averaging, 
the dynamic interaction is momentum independent, and 
thus tC s is momentum independent and separable in the 
energy variables. We therefore write 
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^km,k'm' — T — T K ™ K ™' 

Using Eq. (|), at T = T c the vector 7?, whose km com- 
ponent is K m , satisfies 7? = U dyn * C * 7?, i.e., 



(14) 



The highest temperature for which (jl4|) has a solution is 
T c . Note that Eq. (|l|) is independent of U Btat . Solving 
( |l4| ) numerically, we find 



E 3 
0-32-f , 



(15) 



where we used N = m*/2ir and E F = k 2 F /2m*. Eq. (|l|) 
is in agreement with Ref. Q. As pointed out by Bones- 
teel et al. Q] , T c is reduced by symmetric density fluctua- 
tions, which the composite fermions view as a fluctuating 
pair-breaking magnetic field. Impurities, which exert a 
different disordered potential in each of the two layers, 
further suppress T c . Additional controllable mechanisms 
for T c suppression are discussed below. 

An expansion of the results for small q and u is needed 
for the denominator of JC S in (p"3|), while the numera- 
tor is regular in the long-wavelength low-energy limit at 
T = T c . At finite q and uj m > 0, Eq. (14) is written with 



U dyn and the denominator of C s taken at q and ui m , and 
the sum limited to to' > and to' < —to. The temper- 
ature at which Eq. ( |l4| ) has a solution, which is T c at 
q = u> m = 0, is then expanded for small q and u> (after 
analytic continuation iu> m — > u> + iS). We find, 
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Nq T — T c — iau + r/q 2 
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Numerically we find a « 0.12. The calculation of 77, and 
of the vertex function is difficult since both depend on 
U stat . However, they are related by 



W * J * 7? = 
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(17) 



Here, the expansion of the vertex function to first order 
in q at u> m = = is q^. This equality is a con- 
sequence of a Ward identity. It is shown by writing J 

as 
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(18) 



3 = \ E { d i CS + 2 V * U dyn * V 

fj,=x,y 

+ 2d q ^C s * d q ^U dyD * C s + 2C S * d q ^U dyn * d q ^C 



and comparing with the expansion for 77, described above. 
We note that diagrams with ip^a 2 tp vertices appear in 
( |l8| ) when a derivative with respect to g M is applied to 
a transverse current vertex of U [see Eq. (g)]. Eq. ( |i~7| ) 
leads to a cancelation of both J and 77 from p D . 

Finally, p D may be calculated. The calculation of 
the composite- fermion transconductivity cr* L (|5|), using 
Eqs. ( p6| ) and (0), is similar to the Aslamazov-Larkin 
calculation [[ll] . The conductivity tensor is then inverted 
to obtain the composite-fermion transresistivity, which is 
equal to the (physical) electronic transresistivity p D . As- 
suming p D -C Pxx, we finally obtain Eq. ([!]). We now 
comment on some aspects of that result. 

First, we note that pairing fluctuations lead to an in- 
crease of the transresistivity as the temperature is de- 
creased, in sharp contrast to the temperature dependence 
of the weak coupling contribution. 

Second, the pairing contribution to po is much more 
sensitive to a small density difference An between the 
layers than the weak coupling one. Such a density dif- 
ference amounts to a difference Ap = (dp/dn)An be- 
tween the chemical potentials of the two layers, where 
dp/dn ~ e 2 /ek F . Pairing between the layers is gradu- 
ally suppressed by Ap, in a way similar to the suppres- 
sion of conventional Cooper pairing by spin polarization. 
Complete suppression of the pairing contribution to p D 
happens when A/i ss T c , i.e., when An/n w T C /E F . Sim- 
ilarly, the pairing mechanism is also sensitive to a sym- 
metric deviation of the density in both layers to a filling 
factor v = 5 + Av. In that case both layers are subject to 
a uniform Chern-Simons magnetic field AB oc Av, which 
suppresses pairing. In realistic samples we expect both 
dependences, on symmetric and asymmetric density de- 
viation from v = 1/2, to be mostly determined by the 
inevitable disorder-induced fluctuations in the density. 

Third, we note that the pairing contribution is also 
easily suppressed by the current ji flowing in one of the 
layers. Qualitatively, such a current shifts the composite 
fermion Fermi sphere in the current carrying layer rela- 
tive to the one in the other layer by Ak oc ji, and thus 
suppresses the critical temperature to T c — n(Ak) 2 . In a 
very naive estimate 77 ~ v 2 F jT c and Ak ~ m*j\/ne. 

Experimentally, Lilly et al. Q observed a non- 
monotonic temperature dependence of p D over some nar- 
row ranges of magnetic field in the vicinity oiv = \. The 
precise position of these ranges seems to be sample spe- 
cific, and the non-monotonic temperature dependence is 
observed only at very low currents. The analysis above 
suggests that this phenomenon is the result of pairing 
fluctuations in those parts of the sample at which both 
layers are at filling factor very close to \. The precise 



values of magnetic field at which such parts exist depend 
on the specific disorder configuration. 

A quantitative comparison between theory and exper- 
iment is hard to achieve at this stage, since the one pa- 
rameter in Eq. (|l|), T c , is hard to calculate, and there is 
not enough experimental data to allow for a fit. 

It was recently suggested Q that the low tempera- 
ture behavior of p D observed by Lilly et al. was due to 
the interlayer distance being close to the critical value 
where the two layers form an incompressible state of the 
(1,1,1) type. The mechanism we suggest is different from 
the one leading to a (1, 1, 1) state. In particular, the Hall 
transresistivity is zero in our mechanism, in contrast to 



the finite value predicted in |13|. A measurement which 
separates the longitudinal and Hall components of the 
transresistivity may distinguish between the two mecha- 
nisms. 

In summary, we have calculated the contribution of 
pairing fluctuations to the transresistivity in a system 
of two layers at v = h each. We find that above the 
critical temperature, pairing fluctuations, introduced by 
the interlayer attractive interaction between composite- 
fermions, enhance the transresistivity. In particular, the 
transresistivity does not vanish as the low temperature 
limit is approached. 
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Note added: A related preprint (F. Zhou and Y. B. 
Kim, |cond-mat/980732ip , with which some of our results 
disagree, has been very recently posted on the electronic 
archive. 
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